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Abstract 

In this paper we show that the Schrodinger-Newton equation for spherically sym- 
metric gravitational fields can be derived in a WKB-like expansion in 1/c from 
the Einstein-Klein-Gordon and Einstein-Dirac system. 



1 Introduction 



This work is the sequel to a recent study IU5I1 in which we analysed the qualitative 
and quantitative behaviour of Gaussian wave-packets moving according to the time- 
dependent Schrodinger-Newton equation: 

ihdMt,x)= (-^A-Gm 2 fM^d^^f). (1) 
V 2m J \\x - y\\ J 

In this paper we shall be concerned with the question of how CD can be understood as 
a consequence of known principles and equations. 



We recall that originally CD was suggested in |9l |23|] as a model for the gravi- 
tational localisation of macro objects. This was considered to be an essentially new 
mechanism that also involves hitherto unknown physics. In that sense £D was taken by 
many as a hypothesis to be put to test by experiment. It is then natural to ask whether 
it predicts any observable consequences, say in molecular interferometry, that have 
any chance of being detected in the foreseeable future. It may also be envisaged that 
it might give interesting hints concerning the interface between quantum physics and 
gravitation, and that it might even shed some light on the intricate question concerning 
the necessity of quantum gravity flliril. 



Our results in [15] showed that on the basis of the Schrodinger-Newton equation 
inhibitions of the dispersion of wave packets due to their own gravitational field start to 
become significant at mass scales around 10 10 u for width of 500 nm. This is more than 
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6 orders of magnitude above the current masses in molecular interferometry Ji, 16] but 
only 3 orders beyond the masses already envisaged possible in future experiments. 

This paper is complementary to these attempts and devoted to a better founda- 
tional understanding of the Schrodinger-Newton equation, which we shall view as a 
general model for gravitational self-interaction of matter waves. Accordingly, we shall 
discuss its derivability from known principles, equations, and approximation schemes. 
Through that procedure we also hope to gain insight into the presently justifiable range 
of applicability of £T|), which has been a matter of some debate over recent years. Need- 
less to say that our approach here does of course not exclude the logical possibility that 
CD may capture new physics, as originally anticipated. But even then it would be good 
to know what aspects of it just represent old physics. 

The question we pose closely relates to that of how classical gravitational fields 
couple to "quantum matter", i. e. matter that is described by the Schrodinger equation. 
Note that here "gravitational fields" include those sourced by the quantum matter it- 
self. For classical fields the coupling prescription that derives from the Equivalence 
Principle is the minimal coupling scheme, which proceeds in three steps: 

1 . Formulate the theory of the field you wish to couple to gravity in a way that sat- 
isfies the principles of Special Relativity. In particular: write the field equations 
in a Poincare invariant fashion. 

2. Replace the Minkowski metric by a general Lorentzian metric g and the Levi- 
Civita covariant derivative w. r. t. the Minkowski metric (which is just the partial 
derivative in global inertial coordinates in Minkowski space) by the Levi-Civita 
covariant derivative for the general metric g. 

3. Calculate the energy-momentum tensor T for the matter field from the result of 
step 2 through variation of the action with respect to g and impose on g Einstein's 
gravitational field equation with a right-hand side (energy -momentum tensor) 
including T (and possible external sources). 

There are some ambiguities in step 2 which are well known to result in differences 
in the couplings to the curvature tensor. These are not essential for our concern 
here and will be ignored. In contrast, far more serious for us are the difficulties 
posed by step 1. It requires the formulation of special -relativistic quantum mechan- 
ics, the "non-relativistic limit" (1/c — > 0) of which is Schrodinger's theoryQ How- 
ever, if "special-relativistic quantum mechanics" is inteipreted as "special-relativistic 
one-particle quantum mechanics" this is impossible as such a theory is known not to 
exist (particle creation and annihilation). This seems to force us to interpret "special- 
relativistic quantum mechanics" as "Relativistic Quantum Field Theory". In that case 
step 2 would presumably be achieved for external gravitational fields by considering 
Quantum Field Theory in curved space, but the gravitational self interaction could 
certainly not be captured in that way. 

We conclude that, strictly speaking, there is no obvious way to deduce from the 
standard formulation of the Equivalence Principle the general coupling of a classical 



In this paper we follow standard practise in speaking of "non-relativistic" approximations or limits, 
even though this is unfortunate terminology since it really means to approximate Poincare symmetry 
by Galilei symmetry, both implementing the relativity principle. 
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gravitational field to matter described by Schrodinger's equation. For a spatially ho- 
mogeneous gravitational field g(t) with arbitrary time dependence one could deduce 
the standard coupling, which consists of an addition of V(x, t) = mg(t) ■ x to the po- 
tential in the Hamiltonian, from the standard requirement of equivalence with rigidly 
accelerated frames (see, e. g., ||14|]). This has recently also be shown to work if one 
starts from the Klein-Gordon equation written in co-moving coordinates of a ri gidly 
but time dependent accelerating frame and taking the appropriate l/c — > limit H2QTI . 
But, to stress it once more, the general case does not seem to be deducible by standard 
procedures. 

What we will do is to follow the example just mentioned and formally regard the 
Schrodinger equation as the non-relativistic limit of the Klein-Gordon or Dirac equa- 
tion. To be sure, this is just what is often done in textbooks on Quantum Mechanics, 
and there is indeed nothing to be upset about as long as this is understood strictly as ap- 
proximation scheme for differential equations of complex- valued fields on spacetime. 
However, what remains to be clarified is the interpretation of both fields, and their re- 
lation in the non-relativistic limit. This, we feel, is an important issue that should be 
given due thought before rushing into suggestions to look for possible experimental 
signatures of Schrodinger-Newton equation. In the present paper, however, we con- 
centrate on the formal problem and will only shortly come back to the interpretational 
issues at the end. 



2 The non-relativistic approximation scheme 

We systematise the notion of "non-relativistic approximation" by following the WKB- 
type procedure of Kiefer & Singh 11811 . Recall that the traditional WKB approximation 
is a scheme to obtain the semi-classical limit of a quantum theory by means of a formal 
expansion in terms of a dimensionful parameter H, which is Planck's constant h divided 
by 2-7T. It starts by inserting the ansatz 

?P(x,t) ~ exp Us(x,t)\ (2) 

for the wave function into a given linear partial differential equation, e. g. the 
Schrodinger or Klein-Gordon equation, and subsequently expanding the exponent S 



in terms of H BOH . The equation is then required to be satisfied at each order in h. At 
each order one obtains a certain truncation of the original theory, which makes math- 
ematical sense and may or may not lead to good approximations of the latter. In this 
sense h should be viewed as a deformation parameter of the theory rather than an ap- 
proximation parameter. In the latter case we have to worry about convergence of this 
expansion, where the degree of smallness assigned to h will depend on the context. 
The formal expansion, however, makes sense independent of any context. 



It was shown by Pauli B22I1 how this scheme can be adopted to wave functions with 
multiple components, as for the Dirac equation. Pauli derived as the semi-classical 
limit of the Dirac equation a set of equations that he could not solve in general. This 



was completed by Rubinow and Keller [250 and later by Rafanelli and Schiller [24], 



as well as Pardy 12111 . to yield, e. g., the classical Bargmann-Michel-Telegdi (BMT) 
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equation ^ at the appropriate order in H. An alternative method to obtain the semi- 
classical limit of the Dirac equation using matrix- valued Wigner functions, which also 
leads to the BMT equation, was more recently developed by Spohn J28ll . 

Similarly to the traditional WKB-method, the concept of a non-relativistic limit of 
a relativistic field theory can be understood as appropriate order in an expansion of the 
dimensionful parameter 1/c. Again 1/c should be viewed as deformation parameter 
linking different theories. In particular it allows to contract the Poincare symmetric 
theory to a Galilei symmetric theory as 1/c — > 0. This results in a WKB-like scheme 
applied to the parameter 1/c. Kiefer and Singh [18] showed that this can be used to 
derive the Schrodinger equation from the Klein-Gordon equation, and then generalised 
this method to derive quantum gravity corrections from the functional Wheeler-de Witt 
equation. 

We argue in some detail that this method provides a universal scheme in which h 
as well as 1/c act as deformation parameters, and show that the Schrodinger-Newton 
equation occurs as a 1/c — > limit of the Klein-Gordon and Dirac field coupled to 
Einstein gravity. This goes beyond a former analysis of de Oliveira and Tiomno [Q] 
which uses the traditional method of Foldy-Wouthuysen IU2I1 but does not recouple the 
Dirac field into Einstein's equations. 

Following the strategy just outlined, we can make the ansatz 



(ic^ \ °° 

— S(x,t) J ^2 



(3) 



where tp can be a scalar, vector or spinor field and so are the a n , but S always is a 
scalar function. The semi-classical and non-relativistic limits can then be derived by 
inserting this ansatz into the field equation considered and sorting by powers of either 
h or 1/c. 

We will perform both the semi-classical and the non-relativistic limit for the Klein- 
Gordon equation with an electromagnetic field in paragraph [3] Coupling the Klein- 
Gordon equation to general relativity we show that the Schrodinger-Newton equation 
can be derived as the non-relativistic limit of the coupled Einstein-Klein-Gordon sys- 
tem. 

In paragraph [4] we will repeat this analysis for the Dirac equation with an elec- 
tromagnetic field. We obtain the known results, namely the BMT equation as the 
semi-classical and the Pauli equation as the non-relativistic limit, respectively. When 
coupled to Einstein's equations, the non-relativistic limit of the Dirac equation yields 
the Schrodinger-Newton equation, too. 

Throughout our signature convention for the metric will be "mostly plus", i. e. 
(-,+,+,+). 



3 Klein-Gordon fields 

The free Klein-Gordon equation reads 

' m 2 c 2 \ , 
□ - -nr- 1-0 = 0, 



h 2 



a = -\d 2 + A, 



(4) 



4 



where ip is a scalar field. Introducing an electromagnetic field with electric potential 
and vector potential A by replacing 



d t -> d t + -r4>(x,t) 
n 



dk -> d k - —A k (x,i) 
n 

□ ->■ □ + 



2ie 



fr 2 c 2 
ie 



he 2 



the Klein-Gordon equation takes the form 



(5a) 
(5b) 

(5c) 



2ie 



— f> 2 
h 2 

ie 



2 4 

m c 



+ 



+ c 2 V -A) -0 = 0. (6) 



Note that the last lines of equations (l5cl) and © could be cancelled using the Lorenz 
gauge 4> + c 2 V • A = 0, but we do not want to fix a gauge at this stage because it would 
be c-dependent. 

We now make use of the ansatz ([3]) and calculate the first and second order temporal 
and spatial derivatives for the field ip. They are 



dti> = e 
Vip = e 
d 2 ip = e 
Aip = e 



„2 f 



ic 2 s/h £ ( 



n=0 

„2 00 



ic 2 S//i ^_ ^ [ 
ic 2 s/h c 



n=0 
4 oo 



h 2 

\c 2 s/h ^_ 

h 2 



E 

n=0 
4 oo 



E 

71=0 



Vh 

c 

Vh 

c 

Vh 
Vh 



iSa n + a n -2 
i(V5)a n + Va n _ 2 

— 5 2 a n + 2i5a n ,_ 2 + \Sa n -2 + a n -4 

- (VS) 2 a n + 2i(V5) • Va„_ 2 
+ i(A5)a n _ 2 + Aa„_4 ) , 



(7a) 
(7b) 
(7c) 



(7d) 



where we denote the time derivative dt by a dot and define a„ = for all n < 0. 
Inserting this ansatz into the Klein-Gordon equation © yields 



ic 2 



= exp( T 5 y h 2 



4 oo 



E 



S a„ + 2i5d n _2 + iSa n _ 2 + a n _4 



n=0 \ / 

+c 2 (V5) 2 a n - 2ic 2 (VS') • Va n _ 2 - ic 2 (AS> n _ 2 - c 2 Aa n _ 4 
-^^-Sa n -2 + — - 2e(A • VS 1 )^ + 2ie^4 • Va n _ 2 

e 2 ^ 2 e 2 A 2 2 ie ; • ^ A 

r^-an-4 H ^a n - 2 +m a n + —(pcin-i + ie(V • A)a n _ 2 

h A h h 



(8) 
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3.1 The semi-classical limit 

We first rewrite © by eliminating the h dependence inside the square brackets through 
appropriately shifting the summation index of the terms containing h: 



S 2 a n + 2i5d n _2 + i5 , a n _ 2 + a n -4 



n=0 

+c 2 (V5) 2 a n - 2ic 2 (VS') • Va n _ 2 - ic 2 (A5)a n _ 2 - c 2 Aa n _ 4 

~~5rSa n + -^o n _ 2 - 2e(^4 • V5)o n + 2ie^4 • Va n _ 2 
e 2 <A 2 e 2 A 2 9 ie 



-an H n— On + m a n + ^0a n _ 2 + ie(V • A)a n -2 



i 9 i ^n i o 1 

C C C 



(9) 



Sorting by powers of n we obtain equations 

= ( m 2 - S 2 + c 2 (VS) 2 - - 2e(A • VS) - ^ + ^ 



f J + c 2 V • 1)) a n _ 2 + 2i (5 + J ) 



-2ic 2 [ (VS) - ^2" ] • Va n _ 2 + a n _ 4 - c 2 Aa n _ 4 , (10) 



one for each n. For n = this yields 

= (mc 2 ) 2 - (c 2 S + e(f>y + c 2 (c 2 VS - elY , (1 1) 

which is the Hamilton- Jacobi equation for a relativistic particle. The equations can 
be simplified further if we introduce the four vector 7r u with ttq = —cS — e<p/c and 

7Tfc = —c 2 dkS + eA^. We obtain 

= m 2 c 2 + tt^. (12) 
At order n = 2, now making use of the Lorenz gauge, we get 
= (c 2 S - c 4 AS) a + 2 (c 2 S + e<j^j d - 2c 2 (c 2 VS - eA \ • Va , (13) 
which (with do = dtj c) can be written as 

= (8^)ao + 2^a . (14) 

3.2 The non-relativistic limit 

Again we rewrite equation ©, but this time we eliminate the c dependence inside 
the square brackets through appropriately shifting the summation index of the terms 
containing c. In order not to list terms with index larger than n, we make an overall 
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shift n (n — 2) and compensate for this by an overall multiplication with c 2 /h: 



. 2 \ fi 00 
1C „ \ C° 



c 

n=0 \ 

\2 



S 2 a n -2 + 2i5d n _4 + iSa n _4 + d n _6 



+h(VS) a n — 2ih(VS) ■ Va n _2 — ih(AS)a n _ 2 — 7iAa n _4 
-^^S'a n _4 + ^^d n _ 6 - 2e(^4 • VS , )a n „ 2 + 2iel- Va n _ 4 

H ^— a n -4 + m 2 a n -2 + — (/>a n -6 + ie(V • ^4)a n _ 4 

Sorting by powers of n we now get the equations 

= h{VS) 2 a n 

+ (to 2 - S 2 - iHAS - 2e(A-VS)) a n ^ 2 - 2ih{V S) • Va„_ 2 



(15) 



+— (— ihV — eA) + ( iS -^S ) + 2iSa n _4 



1 2 

--2 (ihd t - e4>) a n _ 6 . (16) 

At order n = this yields simply VS = 0, thus S(x, t) = S(t) depends only on 
time. 

At order n = 2 we then get 

(to 2 - S 2 )a = =>- S = ±mt + const. , (17) 

where the constant term can be ignored and we choose the positive energy solution 
S = —rat. 

Using these results at order n = 4 finally yields the Schrodinger equation 

1 / -A 2 
(ihdt - e<j>) a = — 1-ihV-eA) a . (18) 
2to \ / 

At order n = 6 we get 

(ihdt - &f>) a 2 = — (-mV - el) 2 a 2 - — r {ihd t - ecj)) 2 o . (19) 

Neglecting the vector potential A this reduces to the equation already found by Kiefer 
and Singh flill . Without any electromagnetic potentials, we have 

h 2 h 3 

iha 2 = — - — Aa 2 - - — , AAa . (20) 

2to oto- 5 

3.3 Gravitating Klein- Gordon fields 

Next we consider a Klein-Gordon field coupled to Einstein's equations 

Gp LV = — —T^ u , (21) 



where 



Gfxy — RflV ^ QpuV R 



(22) 



is the Einstein tensor and T^ u the energy-momentum tensor of the Klein-Gordon field, 
the expression of which will be given below (cf. (I38T)). The set of equations © and 
(|2T1) is also known as the Einstein-Klein-Gordon system. We specialise to spherically- 
symmetric metrics which, upon choosing appropriate coordinates, we may write in the 
form H 



(23) 
(24) 



ds 2 = _ e 2A(r,t) c 2 dt 2 + e 2B(r,t) dj ,2 + ^ ^2 + ^2 & ^ 

with determinant 



_ c 2 e 2(A+B) r 4 sin 2 ( 



We expand e A and e B as 



A n (r, t) 



E 

n=0 



An = 1 



Ba = l 



and make use of the further expansions 



-B 



n=0 

oc 



-2A 



n=0 
00 

E 

n=0 

00 



n=0 



C 



G n , 



-2B 



! = E 

n=0 

00 

n=0 

00 

i = E 



n=0 

00 



J B , e - 2B (A'-B') = Y J 



n=0 



Dr, 



Fn 



K„ 



with the coefficients 

Cb = l, 
A> = 0, 
^o = 0, 
^0 = 0, 
Co = 1, 
#o = 0, 
Jo = 0, 

A' = 0, 



Ci 



C 2 = Af - A 2 
D-2 = B 1 — B2 
E2 z 



(A 1 + 5i) ■ Si + B 2 
{A 1 +B 1 )A' 1 + A' 2 



#1 



-2Ai, 



G 2 = 3A? - 2^1 2 

= 3-Bj — 2i?2 



J 2 = -2Ai(ii - fli) - Aiii + B1S1 + i 2 - #2 
K 1 = A! 1 -B' 1 , 



(25) 
(26) 



(27) 
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K 2 = -2J3i(Ai - B[) - AiA[ + B 1 B[ + A' 2 - B' 2 . 



The d'Alembert operator □ in a curved background is 
1 



□ 



-9 



e- 2A c~ 2 ((i - B)d t - d 2 ) + e~ 2B (A' - B')d r 
+ (e- 2 *-l) + A 



(28) 



(29) 



where the dot denotes derivatives with respect to t, the prime denotes derivatives with 
respect to r and A is the flat, three-dimensional Laplace operator. 
The Klein-Gordon-Equation then takes the following form: 



= e- 2A c~ 2 



A-B)%i>-il)) + e~ 2B (A' - B') $ 



+ (e- 2S -l) (V + <) + Atp 



m 2 c 2 



fr 2 



(30) 



Using the same expansion ® for ip as before, and denoting by A r = 2/r d r + d, 
the radial component of the spatial Laplacian, this becomes 



. o \ 4 oo 
1C „\ C 



/ n=0 \ / 



-2A 



A — B) (iSa n _ 4 + d n _ 6 



+S 2 a n -2 ~ 2i5d n _ 4 - i5a re _ 4 - a n _ 6 j + e 2B (A' - B') (i5'a n _ 2 + a' n _ 4 ) 
+ (e~ 2B - 1) (S' 2 a n + 2i5'<_ 2 + i(A r S)a„_ 2 + A r a n _ 4 ) 

Tfl 2 

-(VS) 2 a n + 2i(V5) • Va n _ 2 + i(A,S)a n _2 + Aa n _ 4 ^ a n-2 



(3D 



and with the expansion for the exponentials (|27T ) we obtain 

= exp (jT s ) ¥ £ ( ir) [E{l^(^-^ + 
\ / „=o \ / Lfe=o I 



a n-k-6 



+ ( S a n ^k-2 — 2i5d n _fc_4 — i5a n „. 



(i5'a n _ fc _ 2 + a^_ fc _ 4 ) 



-Hk (—S' 2 a n -k + 2i5 ,/ a^ l _ fc „2 + K^rS)on-k-2 + A r a re _fc_ 4 ) 



m 



-(VS) a n + 2i(V5) • Va n _ 2 + i(AS)a n _2 + Aa n _ 4 — a n _ 2 

n 



(32) 
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As Hq = Jo = Kq = this can still be simplified and for each n we get 

-2 



= h{V ' S) 2 a n + m 2 a n _2 — 2\h(VS) ■ Va ra _2 — i?i(A5)a n _2 — S 2 a n - 



Jk ( i5a n _fc_4 + d n _fc- 



+2i5d n _ 4 + iSa n _4 - fi,Aa n _ 4 + a n _ 6 — 

k=l 

+G k (yS 2 a n _f._2 — 2iS'd n _fc_4 — iSa n _k-4 — a-n-k-6 
+K k {\hS'a n ^ k ^ 2 + ha' n -k-i) 

+H k (-hS' 2 a n _ k + 2\hS' 'a' n _ k _ 2 + ih(A r S)a n - k - 2 + /iA r a n _ fc _ 4 ) 



. (33) 



At lowest order n = (1331 is again equivalent to (VS 1 ) 2 = 0. Thus, S is a function 
of time only and (l33l is trivially fulfilled at order n = 1. 

At order n = 2 we get S" 2 = m 2 and choose, as before, the positive energy solution 
S = —mt. With these results equation (1331 reduces to 

^ A ft 

= -md n _ 4 - — Afl„_4 + — d n _ 6 
2m 2m 

ft 



2~ ^2 [ Gk ( m2a n-k-2 + 2imd n _ fc _4 - d n _ fe _ 6 ) 



fc=i 



+ft (^fea' n _jt_4 + H k A r a n - k -i) - J k {ima n _ k _ A - a n . 



k-6) 



(34) 



At order n = 3 we now obtain G\ = and therefore Ai = 0, G2 = —2A 2 . 
Considering (l33l at order n = 4, we finally see that the Klein-Gordon equation is 
equivalent to the Schrodinger-Newton equation 

h 2 

iMo = — - — Aa + Va (35) 
2m 

with potential 1/ = mhA2- 
3.3.1 Einstein's equations 

Let us now consider Einstein's equations to derive the Poisson equation for the poten- 
tial V. 

The non-vanishing components of the Einstein tensor are 

G tt = e 2A c 2 (\-e~ 2B (\-—Y) (36a) 



1 or 1 2^4' 



G rr = e 2B + + — (36b) 

G tr = — (36c) 

r 

r 2 e -2B( A ,2_ A , B , + A „ + ^ 



e~ 2A \ -B 2 + AB-B) (36d) 
G w = sin 2 0Gg e . (36e) 



10 



From the Lagrangian for the Klein-Gordon field 



C = ~ ({d\){d x r) + IVf ) \/=5 (37) 



2m \ h 2 
the stress-energy-tensor can be derived as 
2 5C 



2m 



m 2 c 2 



(38) 



Its non-vanishing components are 



Tu = ^ <? A m 2 + ^ e^-^'l 2 + ^ (39a) 

T rr = e 2B |</f + I- |V? + ^ e 2 (^ | V| 2 (39b) 

2 2m 2m cr 

Ur = |- + VV) (39c) 
2m 

mc 2 2 I / |2 ^ r2 -2B\ it\2 i ^ r2 — 2j4 i / |2 niiJx 

r» = — g-r M -— e |V I |V>| (39d) 

T w = sm 2 6T e e. (39e) 

We now expand both the Einstein tensor and T^ u /c 4 using Mathematica and 
consider Einstein's equations for each component order by order. We make use of 
the fact that A\ = from our analysis of the Klein-Gordon equation, and we use the 
lower order results to simplify the equations at higher order. The components that are 
not mentioned are trivially fulfilled at the given order. 

n = 0: The it-component yields {rB\)' = 0. 

n — 1: The W-component yields 
2h /3 



.,2 



^-B\ + (rB 3 )'J = 8vrGm|ao| 



(40) 



n = 2: The tt-component yields 

(-4B? - ZrB x B' 2 + (r£ 3 )') = 8vrGm (a> + ofo) . (41) 

The rr- and ##-components both yield B\ = and the ir-component = is then 
trivial as well as the order n = equation. Equations (1401) and (I4TT) then simplify to 

, _ AtrGmr 2 . ,n 

(rB 2 )' = a 2 (42) 

, A-KGmr 2 „ 
(rS 3 ) = (aia + 00%) (43) 
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n = 3: The tt-component yields 

4vrGr 2 



4(A 2 - B 2 ){rB 2 )' + B 2 - 2rB 2 B' 2 + (rB 4 )' 



h 2 



—K\ 2 

2m 1 01 



+ ^ (flo^o - ao«o) + "f/i (^2|ao| 2 + |ai| 2 + a 2 a + a^a 2 ) 



(44) 



The rr-component yields 

B 2 = rA' 2 (45) 

and the 6>#-component is just the derivative of the rr-component. The fr-component 
yields 

-B 2 = 4vriG (a* a' - a a* ') . (46) 

If we define the potential V = mhA 2 as before and analogously U = mhB 2 , and 
also introduce the (r-component of the) probability current 

3kg = y~ {ao a o ~ a*oa' ) , (47) 

we are left with the following set of equations: 

U = rV' (48a) 
(rU)' = 4irGm 2 r 2 \a \ 2 (48b) 
U = -4irGm 2 rj, (48c) 

which are equations for U, V and ao only, together with equations (l43T > and (l44l which 
constrain B%, B 4 and a 2 in terms of U, V and ao- 

3.3.2 Poisson equation 

Let us further analyse the set of equations (|4B1 to show that they are equivalent to the 
Poisson equation for the potential V. 
U is determined from (148bl > to be 

4-irC m 2 l' r 

U(r,t) = — / dff 2 |a (f,t)| 2 . (49) 

Inserting (I48ab into (I48bb yields 



o 



4irGm 2 \a \ 2 = V" + -V = AV (50) 

r 

and therefore the Laplace equation for V that we were looking for. 

We still have to check the consistency of equation (I48cb . Note that by the continuity 
equation which directly follows from the Schrodinger equation we get 

divj+a t |a | 2 = 0, (51) 

where in the spherically symmetric case the divergence is given by 

divj =j' + -j. (52) 

r 
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Differentiating d49l by t then yields 

[/ = 4fGm 2 r f dr T — d t \a (r, t)\ 2 
Jo r 

= — 4-7T G m 2 r dx x 2 divj (xr, t) 
Jo 

= —Gm 2 r \ dV divj(xr, t) 

J unit sphere 

= —Gm 2 r I dSj( r )0 

J 9(unit sphere) 

= -4vr Gm 2 rj. (53) 

Hence, equation (I48cb already follows from the other equations and the system is con- 
sistent. 



4 Dirac fields 

Let us now turn to the Dirac equation and repeat our analysis for this case. The free 
Dirac equation is 

(^-t)*=«- •>•=(! -O-^U' "») <*> 

with the Pauli matrices 

' 1 \ 2 _ / -i \ 3 ( 1 



1 o J' ff =U J' ^=(,0 -1 

V> here is a four component spinor field. 

Denning the matrices /3 = 7 and a k = ^ j k and multiplying equation (1541 by /3c 
yields 

2 \ 

id t + ica k d k - ^-/3 J 1/) = 0. (56) 

Now we can, again, introduce an electromagnetic field with electric potential 4> 
and vector potential A by replacing 

d t ^d t +^(x,t) (57) 
n 

d k -¥d k -'^A k {B,t). (58) 



The Dirac equation then takes the form 



id t + ica k d k -U + ^a k A k - ) ij, = 0. (59) 

an h J 

As for the Klein-Gordon equation, we now make use of our ansatz © with the 
derivatives (Q, where 5 now is a scalar function but the a n are four component spinors. 



2 
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Inserting this ansatz into the Dirac equation (|59l ) yields 



= exp 



lr 2 \ 2 °° 



h h 



n=0 



Sa n + id ra _2 T-a n -2 

n 



-* — > -* TYIC 

-ca ■ (VS)a n + ica • Va n _ 2 + —a • Aa n - 2 — /3a n _2 



(60) 



4.1 The semi-classical limit 

As for the Klein-Gordon equation we rewrite equation (1601 eliminating the h terms: 



■ 9 \ o OO 



5a n + id n _2 itO^ 



n=0 



-ca • (V5)a n + ica • Va n _2 H — a • Ao n — mj3a n 

c 



Sorting by powers of n this yields 



m/3 + S H — g + ca • VS a • A ) a n — id n _2 — ica • Va n _2 = . 



(61) 



(62) 



Making use of the notations ttq = —cS — e(/)/c and tt^ = — c 2 dkS + e^4fc we obtain 
at order n = 



= (mc/3 — 7To — a • tt ) ao 
= (mc - tt m 7^) a 

(mc- 7To)l —a-TT 



ff ■ tt (mc + 7To)l 



(63) 



which has non-trivial solutions if and only if the determinant 

m 2 c 2 — yrg + )(<?• vf) 2 | =0 



(mc — 7To)l —a-TT 

G ■ TT (mc + 7To)l 



vanishes. The Pauli matrices obey the algebra 

a l a j = 5 tj + ie ljk a k (a ■ u)(a ■ v) = u ■ v + iff • (u x v), 



(64) 



(65) 



therefore (a ■ tt) 2 = tt ■ tt and (l64l yields again, as for the Klein-Gordon equation, the 
Hamilton- Jacobi equation for a relativistic particle 



= m 2 c 2 + tt^. 



(66) 



At order n = 2 we get 

(mc - tt^) a 2 = ic7° (d t + ca • a 
= -icV^«o- 



(67) 
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If we name the operators by 

L:=mc- ir^ , D := -icVfy , (68) 

at order n = we have the condition that ao £ Ker(L). At second order we now have 
Dao G Im(L) which is equivalent to Dclq € (Ker(L' )) . Now note that = 7°L7° 
and therefore x S Ker(L^) 44> 7°x S Ker(L). 

The condition that Da^ is in the image of L is therefore equivalent to the condition 
that for any two solutions ao, So £ Ker(L) to the first order equation we have 

aol^d^ao = 0, (69) 

where ao = (ao) ' 70 is the adjoint spinor. Equations (l63l and (l69l together determine 
the solutions ao at first order. 



4. 1. 1 Derivation of the BMT equation 

We can use these results to obtain the Bargmann-Michel-Telegdi equation |Q] as a 
necessary condition. First multiply equation (I67T ) by (rnc + ^7^) from the left. Then 
the left-hand side vanishes and we get 

(mc + n^)Yd„a = 0. (70) 

Now let us calculate 

-2vr^a = ^(-2^)^ao 

= -mcrfd v aQ + d v (7 !/ ^7 M «o) - (d^^^ao 



(d^ L - a ■ Q# - Vvr ^ + i7 5 a '(vxjfjj 
((9„7r M )a - e ( -a ■ ( E - A) + i>y 5 a • B ) a 



ao 



c 

= ( V> + fa-f-fao, (71) 

where 7 s = ij ^ 1 ^ 2 ^ 3 and we made use of (PTOl in the third and (1631 in the fourth 
line. Performing the same calculation for the adjoint leads to 

- 27r^ao = (d^)a - | F^ao-f-f . (72) 



We follow t2M and set 

S" = a o7 Vao, (73) 
the spin density, then we get from equations (1711 ) and (1721 

tt^S" = ^(9 M ao)7 5 7"a + a o7 5 7^^ M a 

= -(d^ST + ^a ( 7 p 7 CT 7V - 7 5 7*W) «o 

= -(d^)S u + \F pu a^ (2^ 7 CT - 2^ 7 p ) a 

= -(3^)5" + eF%5". (74) 
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For the normalised spin density 

s * - OT S " (75) 

the first term vanishes and we obtain the BMT equation (for g-factor g = 2) 



d - 1 e 
—S u = —^daS v = - 
dr m m 



S u = —^d^S v = —F v p S p , (76) 



where the first equality holds because 7r M = mv^ with the four-velocity of the 
relativistic particle in an electromagnetic field. 

4.2 The non-relativistic limit 

Again, we rewrite equation (1601 eliminating the c terms: 



2 \ 3 00 

ic „\ cr 



= exp (T s J^E( v 

' n=0 \ / 



+ lOn-3 ^ a n-3 

ft 



(77) 



—Vha ■ (VS)a n + iv^a • Va„_2 H — pa • Aa n -2 — m(3a n -i 
Sorting by powers of n we get 

Vha ■ (VS)a n + (jS + a n _i + ~y=® ' i^V — eAj a n _2 

- - (ifi^ - &j>) a n _ 3 = 0. (78) 

At order n = this simply becomes VS* = 0, and therefore S = S(t) is a function 
of time only. 

Now, splitting the four component spinors a n = (a Hi i, a n> 2, a n ,3, a n ,4) i nto two 
two component spinors a> = (a n ,i,a n ,2) and a< = {a n ^,a n ^), at order n = 1 we 
get the two equations 

(m + S)a£ = (79a) 
(m - S)a£ = 0, (79b) 

which can be consistently fulfilled only if either S = —mt and the negative energy 
component vanishes or S = +mt and the positive energy component vanishes. 
From here on, we will choose the first and set Oq = 0. With this choice, the Dirac 
equation (1781 simplifies to the following two equations: 



= Vha ■ (-iW - eAj a<„ 2 ~ (V&t ~ a>_ 3 (80a) 

a «-i = o 1 m d ■ (-^V - el) a>_ 2 - — - e0) a<„ 3 . (80b) 

At order n = 2 (I80ab is trivially fulfilled. Equation (I80bb determines to be 

af = — l —=a-(-\hV-eA)a^. (81) 
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Using this at order n = 3 in equation (I80ab . we get 



1 

2m 



2 



a ■ -ihV - eA 
\2m V \ 

Making use of 

f-iW - eAj x f-iW - eA\ = -h 2 V x (Va^) 



(ihd t - e0) Uo = °- ( 82 ) 



=o 

> 



+ ieh V x (a5*A) +\ehA x Va^ + e Ax A = iehBafi (83) 



=ag'VxA- J 4xVa5' 



and the algebra (l65l l of the Pauli matrices we end up with the Pauli equation 

1 / -a 2 eh ^ 

(ihd t - e6) = — -iW - a^" - —a • Ba^. (84) 



2m V 7 u 2m 

Equation (1 80b I) at order n = 3 determines ajf to be 

= - — a ■ (-ihV - eAj . (85) 

This is exactly the same relation as (f8TT > with the indices shifted by one. Thus, from 
equation (180al) at order n = 4 we will get the Pauli equation again for a>. But equation 



80b I ) at order n = 4 now has an additional term depending on af 

a ■ f-ihV - el) a% + — (-ihd t + e<p) af . (86) 



2mVh V / 2m 

This will make the equation for more complicated, and different from the Pauli 
equation, at order n = 5. 

4.3 Gravitating Dirac fields 

Now let us consider the Einstein-Dirac system, consisting of the Dirac equation (1541 
and Einstein's equations (|2~TT) . 

As for the Einstein-Klein-Gordon system, we make an ansatz using a spherically 
symmetric metric tensor and therefore also need the stress- energy-tensor to be spher- 
ically symmetric. Note that a single Dirac particle cannot have spherical symmetry, 
which is why we have to average over all spin directions. But we will only take this 
into account at the very end of our considerations. 

We use the general relativistic formulation of the Dirac equation according to Fin- 



ster niQTI and Finster et al. Ill ill 



(p - ™) V = 0. (87) 

The Dirac operator is defined as 

Jj) = iV(x)d fl + Y(x) (88a) 

Y(x) = T^x)Z^x) (88b) 

Z„ = \pd^p - ^Tr (r"V M r') r„r p + ^Tr (pr M V,n p (88c) 

P = ^ upa T^T^, (88d) 
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where V is the covariant derivative, e pupcr is the Levi-Civita symbol defined by 

+1 if (pvpa) is an even permutation of (trOip) 
t-ixvpa = yf—g \ — 1 if (fiupa) is an odd permutation of (tr9ip) (89) 
if two or more indices are equal 

and T M is a representation of the space-time dependent Dirac matrices, satisfying the 
Clifford algebra 

{T^V} = -2g pu . (90) 

The additional part Y(x) will turn out not to contribute to the Schrodinger-Newton 
equation and is only relevant at higher order in 1/c. 

It is useful to represent these matrices in the basis where they become the linear 
combination 

r* = e -A c -1 7° (91) 
r r = e~ B (7 1 cos 9 + 7 2 sin 9 cos 9? + 7 3 sin 6 sin (p) (92) 

T =-(—7 sin 9 + 7 cos 9 cos p + 7 cos 9 sin pj (93) 



- f-7 1 sin# + ~ 2 
r 

= — !— (-7 2 sin 99 + 7 3 cos 99) , (94) 

of the Dirac matrices 7^ as defined before. These matrices satisfy the anti-commutator 
algebra, and simplify the equations because in this representation p = 7 s and there- 
fore the first term of Z p vanishes because p is constant, and the third term vanishes 
because derivatives of the Y p as well as the T p themselves are linear in the 7^ and 

Tr ( 7 5 7^7 V ) = 0. Therefore 

y = - y^Tr (r*v M r") t»t v t p (95) 
= -^Tr (r"v„r>) (^r„ - 6£r p - v^ 9vp + leW) (96) 

= _ i T r (TVpT") + ie^Tr (r„ V^,) r CT7 5 (97) 

= -^(r%r")r„ (98) 

o 

All this is in agreement with JlHl and can be straightforwardly verified. The third line 
follows, because the by terms in the second line are equal and the g up term vanishes. 
The second term in the third line vanishes, because the trace vanishes if all three indices 
are different. 

Using that V p T p = a p T p is some linear combination of the gamma matrices we 
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get 



TV (T u V p T p ) = a p Tr (T V T P ) 



and therefore 



a p Tr({T»,T p } + [r\T p ]) 
^a p (Tr ({T u , F p }) + Tr (FT") - Tr (FT")) 

^a p (Tr (-2^1 4x4 ) + Tr (FF) - Tr (FT')) 

-a p ^Tr(l 4X 4) 
-4a" 



= iF3 M + -v M F. 



For our ansatz for the Dirac matrices we get 
1 



v M F 



-9 



1 



i 



=<9,. 



Br* 



+ 



1T») + 



ia M f. 



coter*-2-r r 



7^ cos + sin (g 
r sin 



BT t +[A' + -)r r -—r 



r sin l 



(7 1 cos 9 sin — (1 + cos 2 9)(^f 2 cos ip + 7 3 sin</?)) 



sF + (A' + - (1 - e s ) r 



and therefore 



B 



Jt> = irM d f + - + if [d r + — + 



A' 1 - e 



B 



+ iFa, + iFd ( 



(99) 
(100) 

(101) 

(102) 

(103) 
(104) 

(105) 
(106) 



(107) 
(108) 

(109) 



Now we make use of the expansion Q for if) as before and insert it into the Dirac 
equation (I87T ). We introduce the matrices 5 analogously to the previously defined a. 



a r = 7°7 1 cos 6 + 7°7 2 sin 9 cos ip + 7°7 3 sin 9 sin ip 
a e = — 7°7 1 sin 9 + 7°7 2 cos 9 cos <p + 7°7 3 cos # sin ip 
= — 7°7 2 sin 93 + 7°7 3 cos <p, 



as well as the rotated Pauli matrices 



a 



(7 



a 1 cos 1 



a 2 sin 9 cos 99 + <t 3 sin 9 sin (/? 



-o" 1 sin 9 + a 2 cos cos ip + a 3 cos sin ip 
-a 2 sin 93 + <t 3 cos 



(110a) 
(110b) 
(110c) 

(111a) 
(Hlb) 
(111c) 
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Multiplying the Dirac equation (1871 ) with 7 = /? from the left we then get 

■ /— \ '• 



• 9 \ O OO 

1C „\ c 



= exp| T s ; * 



n=0 



1 _ 4 a i _ /i . 

e ^a n _i + — e a n _ 3 



l) a r S'a n + ia • Va n _ 2 



+i (e" B - 1) a r a' n _ 2 + ^e" jB A'a r a n _ 2 + - {e~ B - l) a r a n ,_ 
m 

and with the expansion for the exponentials (1271 ) we have for each n 

= v^a • (VS)a n + /3ma n ,_i — ix/^S • Va n _2 



(112) 



+£ 

fc=0 



C k Sa n -k-x - iC k a n - k -3 - -E k a n _k-3 + VhD k a r S'a n _ k 



-iy/hD k a r a' n _ k _ 2 



——l< k a a n - k - 2 D k a a n -fc-2 



(113) 



V^a • iyS)a n + f3ma n -i + Sa^-i — iv^a • Va n - 2 



+£ 



fe=i 



VhD k a r S'a n _ k + C fc 5a n _ fc _i - iVhD k a r a' n _ k _ 2 



Fk + -D k ) c/a„-fc_2 — iCfca n _fc_3 — -E k a n _ k -z 
r 2 



(114) 



As before, at order n = this yields V5 = 0, i. e. S is a function of time only, and 
at order n = 1 we get (/3m + 5) an = 0. Again, we split the four component spinors 
o-n = («n,l> On,2 5 «n,3> a>n,4) into two two component spinors a> = (a n; i,a n) 2) and 
a< = (a n> 3, a TO) 4) and choose Oq = 0. We are then left with S = —mt and from (II 141 ) 
we get the following two equations: 



n 

= iVha ■ Va<_ 2 + id>_ 3 + ^ ^mC k a^ i _ k _ 1 + \\fhD k a r a, 



< 1 

n-k-2 



k=l 



"n-1 



+ 



2 
2m 

4m 



(115) 



a ■ Vo 



i . < <A r 1 n < , iv^ n ~ r > / 
2 " 2^ a ™- 3 " A. [a * "-*-! + 2m fcff a "- fe - 2 



fc=i 



F fc + -D k ) a r a>_ k _ 2 + — C k a<_ k _ 3 + — £ fe a< 



4m 



jfc-3 



• (116) 



At order n = 2 the first equation yields = mCia^ and therefore C\ = 0. This 
means A\ = 0, C 2 = —^2 and i*\ = 0. The second equation yields 



iVK 
2m 



cr • Van • 



(117) 
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Inserting this into the first equation at order n = 3 we get 



= h/TSr- V ( Va5* j +ia£ -mA^ (118) 



U = W Iw ■ V I — 

As (of • u) 2 = for any vector u we obtain the Schrodinger-Newton equation 



h 2 

iM> = -—Aa> + Va>, (119) 

with the potential V = mhA-i as we did for the Klein-Gordon equation. 
The second equation at order n = 3 yields 

< = _^.Vaf + ^(V + ^). (,20, 

This result could now again be inserted into the first equation (11151 ) at order n = 4 
which would result in the evolution equation for a^, but in contrast to the pure Dirac 
equation (section l4~2l where we obtained the Pauli equation also for af, this evolution 
equation will be different from the Schrodinger-Newton equation. 

4.3.1 Einstein's equations 

Let us first derive the stress-energy-tensor for the Dirac field. The Dirac-Lagrangian is 

/ Trie \ 

C = hc4[Ip - —)^yf^. (121) 



h ) 

Considering the variation with respect to g^ v on-shell (i. e. assuming the Dirac equa- 
tion to be satisfied) we have 

5C = he Re \$ (i (oT^) + SY) ^V^d] • (122) 



According to Finster et al. Illlll . with our special choice (1911 for the Dirac matrices 



Re (#YV0 = \^ VIXJ ($9ux) Tr (r A ^T p ) (123) 

5T» = ~g^ (5g uX )T x (124) 

Sg^ = -g^\gv^g Xa (125) 
and therefore the stress-energy-tensor is 

he 

r M „ = -— Re [^(il^ + xY^) xl>] 

he — 
+ v e b^Tr (T u d a r p ) + g u \Tic (T„d a T p )] ^T^. (126) 
o 

Note that (11241 ) is determined up to local Lorentz transformation. 5T^ has to obey 

8{T' i ,T v } = -6g tu ' => {ST' i ,r v } = -Sg^ (127) 
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and (11241 ) is a special solution to this equation. 

The non-vanishing components of the stress-energy-tensor for a Dirac field are 
then given by 



T 

Ttr 

Tte 



-hcRe -ie 2A c 2 V>r*^ 

-HcRe [ie 2S ^r r V] 
He 



2 

he. 



Re 



-Re 



2 

he* 

2 



-ie 2A c 2 #V + ie 2B i>T r ip 
ir 2 ^T e i) - ie 2 V^r*<%V 



J- t^p 

T r e 
T 

- 1 rip 

Tee 
T 



he 



ir 2 sin 2 diffVip - ie^c 2 ^*^ 



hf 

2 



l-e B \- 



hc 

y Re 



ir 2 v>rV + ie 2B ^r r ^v 



h 

"2 6 



(128a) 
(128b) 
(128c) 

(128d) 
(128e) 



B-AJl 



r 1 sin 9 BifyfTfif) (128f) 



.^Re [ir 2 sin 2 fl^ty' + ie 25 ^^] + |e B " A r 2 sin OBijyfT V 



-/icRe 



ir 2 ^r e a^ 

ir^^d^ip + ir 2 sin 2 9^T v dgil> 
-hcRe [ir 2 sin 2 fl^r^V] • 



- y Re 



Up to order c° these can be simplified to 



T tt = 
T — 

± rr — 
Ttr = 

T te = 
T r e = 
Tt v = 
T — 



-hcRe 

-hcRe 
He 

- y Re 

he 
"T 

-^Re 
2 

he 

—Re 

he 

"Y 



V' 1 ' (-ie A c) </> 
(ie B 5 r ) V' 
(-ie A c) V' + V^ (ie B 5 r )^ 



( vroP ) i/> 



To 



= —hcRe 

= -^Re 
2 



T w = -hcRe 



+ O(c ) 

$ (\ra e ) %// + $ (ie B 5 r ) dgip +O(c ) 
^ (ir sin 6a v ) ip] + O(c ) 

(ir SV0.6W) if>' + V f (ie B 5 r ) + 0(c°) 

ir sin 9a v ) dgip 

ijj' (ir sin Oa^) d^ip . 



(128g) 
(128h) 

(128i) 
(128j) 

(129a) 
(129b) 

(129c) 
(129d) 
(129e) 
(129f) 
(129g) 
(129h) 
(129i) 
(129j) 
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Again, as for the Klein-Gordon stress-energy tensor, we expand Tj^/c 4 using 
Mathematica and consider Einstein's equations with the Einstein tensor given in equa- 
tion (1361) for each component order by order, making use of A\ = 0. The result differs 
from the Klein-Gordon result only slightly. 



n = 0: The it-component yields (rB\)' = 0. 



n — 1: The ii-component yields 

n = 2: The M-component yields 

2fi 



>|2 



^ (-4B? - 3rBxB' 2 + (rS 3 )') = 8nGm (of f a5" + a^ T af ) . 



•>t„> 



(130) 



(131) 



The rr-components yields Si = 0, the ir-component B\ = and the ##-component 
B[ = are then trivial as well as the order n = equation. Equations (11301 ) and (11311) 
then simplify to 

AirGmr 2 , 



(rB 2 )' 



,>|2 



, 47rGmr 2 / >t > >t > 
(T-°3j — 1 l a l a + a n a 



A "1 



(132) 
(133) 



n = 3: The M-component yields 



4(^ 2 - B 2 )(rB 2 )' + Bl- 2rB 2 B' 2 + (rB 4 )' 



4vrGr 2 

— w 



ih ( >t-> 



"0 a 



+ mft ^yl2|aol 2 + l a i| 2 + a 2^ a o > + a ^ a 2j 
B 2 = rA' 



(134) 
(135) 



The rr-component yields 
and the ##-component is just the derivative of the rr-component. The fr-component 

8vrGm, 



yields (with a r as defined in (11111) ) 

-B 2 = 2n\G (a^a^' - a^'^a^ 

'„>t„>' ,>'L>) , A T OT ,>t 



. Re ( a>Q^a r af 



2tt'xG ( Oq Oq — a,Q 'Oq J + 47rGRe iV + axV 



4vriG ( a^ j a>' - a^'^a> ) + 4vrGRe 



>'t„> 



(136) 



where we made use of (11171 ) to express af in terms of Oq in the second line. The t6- 
and ^-components yield 



= Re (apa 9 af J <^ = Refa^ 1 " (iV + a x VJ^ ajf] (137a) 
= Re (apa^af ) O = Refa^ 1 " (iV + a x v) a^] . (137b) 
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Now, because we consider the spherically symmetric Einstein's equations we have 
to average over all spin orientations in order to get a symmetric stress-energy-tensor. 
This implies that terms proportional to a x V in equations (11361 ) and (11371 ) will van- 
ish and that (11371 ) is equivalent to vanishing 6- and ^-components of the probability 
current 

Joirac = ^ ((Va>t) a> - a>^Va>) . (138) 

If then, again, we define the potentials V = mhA2 and U = mhBi and define j 
as the r-component of the probability current we get the same set of equations as in 

(EH) 

U = rV' (139) 

(rU)' = 4vrGmV|a5 > | 2 (140) 

U = -AirGm 2 rj, (141) 

together with equations (11331 ) and (I134I ). As shown in subsection I3.3.2l this yields the 
Poisson equation for the potential V. 



5 Summary 

From our analysis of the Klein-Gordon and Dirac equations we conclude that the ex- 
pansion in either H or 1/c according to our ansatz © for the fields is a valid scheme 
to obtain both semi-classical and non-relativistic limits of field equations in an unam- 
biguous way. Applying the same ansatz to the self-gravitating Klein-Gordon and Dirac 
fields, as mathematically represented by the Einstein- Klein-Gordon and Einstein-Dirac 
systems, leads to the Schrodinger-Newton equation. Hence we may say that the Schro- 
dinger-Newton equation follows from the self-gravitating fields in the same way as the 
linear Schrodinger equation can be derived in flat space. Seen from that direction one 
concludes that the Schrodinger-Newton equation should provide a better description 
than the linear equation. 

However, one may ask: description of what? In order to arrive at our result we 
considered the classical (i. e. not second quantised) Klein-Gordon or Dirac field as 
source for the classical gravitational field. The Schrodinger function then merely ap- 
pears as part of these fields in the 1/c expansion and with certain phases (due to the 
rest mass) subtracted. The central question is, whether this is the right way to represent 
the gravitational field of quantum systems. Reading the classical fields as one-particle 
amplitudes, it amounts to assuming the validity of the semi-classical Einstein equations 

Rp,v - \Rg^ = (142) 

where now represents the full (second quantised) field. This equation has a long and 
controversial history, which we will not review here. The interested reader is referred 
to the book by Kiefer fl7ll a nd references therein, and the recent discussion in Carlip 
fl, Salzman and Carlip ibfill . Interesting early (1957) discussions between Feynman 
and others on that subject may be found in [H]. 

An obvious objections is this: If there were a gravitational self-coupling of the 
Schrodinger function, then there should also be an electromagnetic self-coupling. In 
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a lowest 1/c expansion (neglecting magnetic fields) this would again result in equa- 
tion CD) with Gm? replaced by — e 2 /(4-7reo), leading to enhanced dispersion due to 
electrostatic repulsion. Applied to the hydrogen atom, where the electron now not only 
"sees" the electrostatic attraction of the proton but also the electrostatic repulsion of 
its own charge-cloud, it seems obvious that this cannot again reproduce the known en- 
ergy spectrum. Interestingly, precisely this idea of implementing the electromagnetic 
self-interaction of the quantum-mechanical wave function occurred to Schrodinger im- 
mediately after he wrote his famous papers on wave mechanics. In 1927 he argued J27I1 
that, as a matter of principle, the self-coupling was required by consistency in order 
to get closed systems of field equations. Just like one derives the radiation reaction 
of charges through interaction with their own field in ordinary Maxwell theory. But 
being convinced that the result of this was incompatible with observed facts, like the 
hydrogen energy levels, Schrodinger concludes that, from a classical field-theoretic 
point of view, there is a "strange violation of the closedness of the field equations'!!. 

Schrodinger's idea has been revived in the mid 1980's by Barut and collaborators, 
who wrote down a rather obvious non-linear Dirac equation which is obtained from 
the linear equation by first splitting the field into a self- and external part and then 
eliminating the self-part by means of Green's functions and the self-current, the latter 
then introducing non-linearities. They claimed that without any further input from 
QED this suffices to account for effects usually attributed to the quantisation of the 
electromagnetic field, like Lamb shift, spontaneous emission, and anomalous g — 2. 
And even though from a path-integral perspective their equation is, in fact, just that 
one obtains from an effective action after integrating out the photons, it seems very 
far reaching indeed. What is more, it was also suggested (but not shown) that the 
correct values of the hydrogen spectrum can be obtained if due care is taken of the 
non-perturbative nature of this self interaction 

Another problem concerns the application of the Schrodinger-Newton equation to 
molecular interferometry. What aspect of the complex molecule, comprising many de- 
grees of freedom, does our tp represent? It just depends on 3 coordinates, so one might 
be tempted to identify it with the centre-of-mass motion. But one would certainly not 
expect an influence of the centre-of-mass motion from gravitational pair-interactions 
of many -particle systems yj]. Self-interactions are also absent in a manifest Galilei 
invariant (i.e. manifest "non-relativistic") second-quantised theory of the Schrodinger 
field in a Newton-Cartan spacetime |0]. This still leaves open the possibility to think 
of the Schrodinger-Newton equation as fundamental deviation from one-particle quan- 
tum mechanics. Alternatively one can regard our if) as a condensate and think of the 
non-linearity as coming about through effective interactions, just like in the Choquard 



equation IU3L I19H. 

Our results concerning the derivability of the Schrodinger-Newton equation by 
WKB-like methods from the Einstein- Klein-Gordon and Einstein-Dirac systems are, 
as such, independent of these questions. But, as stressed in the introduction, they 
clearly need to be addressed, in one form or another, in any attempt to make a well 
founded physical applications. 



2 "Gerade die Geschlossenheit der Feldgleichungen erscheint somit in eigenartiger Weise durchbrochen" 
(p. 271 (H). 
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